function sys=plant(t,x,dt,Ft)
sys=zeros(2,1);

x1=x(1);x2=x(2);eps=1;

sys(1)=eps*(x1-(x1^3)/3+dt)+x2;
sys(2)=-x1+Ft;
end